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Abstract 



We study the resummation of large logarithmic perturbative corrections to the partonic 
cross sections relevant for the process pp —>■ hX at high transverse momentum of the 
hadron h. These corrections arise near the threshold for the partonic reaction and are 
associated with soft-gluon emission. We perform the resummation to next-to-leading 
logarithmic accuracy. We present numerical results for the fixed-target regime and 
find enhancements over the next-to-leading order cross section, which significantly 
improve the agreement between theoretical predictions and data. We also apply the 
resummation for RHIC kinematics and find that subleading terms appear to play a 
rather important role here. 



1 Introduction 



Cross sections for single-inclusive hadron production in hadronic collisions, H1H2 hX, play an 
important role in QCD. At sufficiently large hadron transverse momentum, pt, one expects that 
QCD perturbation theory can be used to derive predictions for the reaction. Since high pt implies 
large momentum-transfer, the cross section may be factorized at leading power in pt into convolu- 
tions of long-distance pieces representing the structure of the initial hadrons and the fragmentation 
of a final-state quark or gluon into the observed hadron, and parts that are short- distance and 
describe the hard interactions of the partons. The long-distance contributions are universal, i.e., 
they are the same in any inelastic reaction, whereas the short-distance pieces depend only on 
the large scales related to the large momentum transfer in the overall reaction and, therefore, 
can be evaluated using QCD perturbation theory. Because of this, and because of the fact that 
single-inclusive hadrons (e.g., pions) are rather straightforward observables in experiment, cross 
sections for H1H2 — > hX offer a variety of important insights into strong interaction dynamics. 

If the long-distance pieces, parton distribution functions and fragmentation functions, are 
known from other processes, especially deeply-inelastic scattering and hadron production in e^e~ 
annihilation, one may test the perturbative framework outlined above. In particular, one may 
examine the relevance of higher orders in the perturbative expansion. Any discrepancies between 
the predictions and experimental data may also provide information about power-suppressed con- 
tributions to the cross section. 

Alternatively, one may also gain information about fragmentation functions. For example, 
e~^e~ annihilation is mostly sensitive to quark-to-hadron fragmentation functions, whereas data 
from hadronic collisions may also provide information on gluon fragmentation. In addition, the 
reaction H1H2 — »■ hX may be used to probe the structure of the initial hadrons. Of particular 
relevance here are spin effects, associated with polarized initial protons. At the BNL Relativistic 
Heavy-Ion Collider (RHIC), one measures spin asymmetries in polarized pp hX scattering, in 
order to investigate the spin structure of the nucleon. Finally, high-p^- hadrons are also important 
probes of strongly interacting matter in a high-energy nuclear environment, as generated by heavy- 
ion collisions. Here, hadron production in proton-proton collisions provides an important baseline 
for the study of nuclear dynamics. 

Whatever the uses of processes H1H2 hX, a central piece is in each case the perturbative 
partonic hard scattering and our ability to reliably evaluate it. Lowest-order (LO) calculations of 
the partonic short-distance cross sections were performed a long time ago pQ , and later improved 
when the next-to-leading order (NLO) corrections were computed On the experimental 

side, an extensive data set on high-p-r single-inclusive hadron data has been collected, both from 
scattering off fixed targets and from colliders at much higher energies. Detailed comparisons 
of NLO calculations with the experimental data have been carried out recently in [3 13 El- 
They show the overall trend that NLO theory significantly underpredicts the cross section data 
at fixed-target energies, but yields a good description [HI El HD] of the collider data. 

In the present paper, we further improve the theoretical calculations by implementing the 
all-order resummation of large logarithmic corrections to the partonic cross sections. At partonic 
threshold, when the initial partons have just enough energy to produce a high-transverse momen- 
tum parton (which subsequently fragments into the observed hadron) and a massless recoiling 
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jet, the phase space available for gluon bremsstrahlung vanishes, resulting in large logarithmic 
corrections to the partonic cross section. To be more specific, if we consider the cross section 
as a function of the hadron transverse momentum px, integrated over all hadron rapidity, the 
partonic threshold is reached when \/§ = 2pT, where a/S is the partonic center-of-mass (cm.) 
energy, and pT = pr/z is the transverse momentum of the produced parton fragmenting into the 
hadron, the latter taking the fraction z of the parton momentum. Defining Xt = 2pt/\/§, the 
leading large contributions near threshold arise as In^'^ (1 — x^) at the kth order in perturbation 
theory, where as is the strong coupling. Sufficiently close to threshold, the perturbative series will 
be only useful if such terms are taken into account to all orders in ag, which is what is achieved 
by threshold resummation ^1] ^1 . This resummation has been derived for a number of cases 
of interest, to next-to-leading logarithmic (NLL) order. As far as processes with similar kine- 
matics are concerned, it has been investigated for high-py prompt photon production in hadronic 
collisions [HJ UHl HEl E], and also for jet production [T21 HH CH], which proceeds through the 
same partonic channels as inclusive-hadron production. We will actually make use of the results 
of |18[ for the resummed jet cross section in our analysis below. 

The larger xt, the more dominant the threshold logarithms will be. Since s = X1X2S, where 
Xi^2 are the partonic momentum fractions and is the hadronic cm. energy, and since the 
parton distribution functions fall rapidly with increasing xi^2, threshold effects become more and 
more relevant as the hadronic scaling variable xt = 2pT / \/S goes to one. This means that the 
fixed-target regime with 3 GeV ^ Pt ^ 10 GeV and of 20—30 GeV is the place where threshold 
resummations are expected to be particularly relevant and useful. We will indeed confirm this 
in our study. Nonetheless, because of the convoluted form of the partonic cross sections and the 
parton distributions and fragmentation functions (see below), the threshold regime xt — ^ 1 plays 
an important role also at much higher (collider) energies. Here one may, however, also have to 
pay attention to terms that are subleading near threshold. 

In Sec. 121 we provide the basic formulas for the inclusive-hadron cross section at fixed order in 
perturbation theory, and display the role of the threshold region. Section IHl presents details of the 
threshold resummation for the inclusive-hadron cross section. In Sec. ID we give phenomenological 
results. We focus primarily on the fixed-target regime, but also give some exploratory results 
for collider energies. We do not present an exhaustive phenomenological analysis of all hadron 
production data available, but select some representative examples. Finally, we summarize our 
results in Sec. The Appendix compiles some useful formulas for the threshold-resummed cross 
section. 



2 Perturbative cross section and the threshold region 

We consider single-inclusive hadron production in hadronic collisions, 

Hi{P,) + H2{P2) h{Ps) + X , (1) 

at large transverse momentum pt of hadron h. We integrate over all angles (equivalently, pseu- 
dorapidities rj) of the produced hadron. We note from the outset that this does not directly 
correspond to the experimental situation where always only a certain range in pseudorapidity is 
covered; we will return to this point later on. The factorized cross section for the process can then 
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be written in terms of the convolution 
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where f/ is the pseudorapidity at parton level, with 57+ = — r)_ = In 

sum in Eq. runs over all partonic subprocesses ab —>■ cX, with partonic cross sections ddah^cx-, 
parton distribution functions fa/Hi and fb/H2i aiid parton-to-hadron fragmentation functions Dh/c- 
The scales iipi and ^ff denote the factorization scales for the initial and final states, respectively. 
The dependence on them, and on the renormalization scale /ir, is implicit in the partonic cross 
section in Eq. (j2I). 

The partonic cross sections are computed in QCD perturbation theory. Their expansions begin 
at 0{a\) since the LO partonic processes are the 2 — 2 reactions ah cd. Therefore, 



daab^cxixl, fj) = alifiR) rfo-f^Ldl^r, v) + «s(/iR) da^aLcxi^T, v) + C(«s) 
It is customary to express x^ and fj in terms of a different set of variables, v and w: 



Xrp = 4:vw{l — v) 
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At LO, one then has 
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where the 5(1 — w) function simply expresses the fact that xrcosh(i7) = 1 for 2 — > 2 kinematics. 
It allows to trivially perform the 7) integration of the partonic cross section. Defining 



V ^-2^- r+ x^sdaab^cx{xl,ri) 
the LO cross section for the process gg gg becomes, for example. 



(6) 



y(o) /^2n 



~2 N 3 

1 ±r_ 

4 
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Analytical expressions for the NLO corrections c?a 
ically, they read: 
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have been obtained in Schemat- 
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where the + distributions are defined in the usual way. 
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The function F{w,v) in Eq. (jSJ represents all remaining terms without distributions in w. The 
terms with + distributions in Eq.® generate at NLO level the large logarithmic contributions 
we discussed earlier, and which we will resum to all orders in as- After integration over 17, 
the term [ln(l — w)/{l — yields a contribution oc ln^(l — to ^ab^cx^ P^^^ terms less 
singular at xt = 1- At higher orders, the leading logarithmic contributions are enhanced by terms 
proportional to [ln'^''~^{l — w)/{l — w)] + in da^l^^xi'^^y w)/dv dw, or to a| ln^*^(l — x^) in S^^|_^^j(.. 
As we discussed earlier, these logarithmic terms are due to soft-gluon radiation and, because there 
are two additional powers of the logarithm for each new order in perturbation theory, may spoil 
the perturbative expansion unless they are resummed to all orders. 

As follows from Eq. (j21), since the hadronic variable xt is fixed, xt assumes particularly large 
values when the partonic momentum fractions approach the lower ends of their ranges. Since 
the parton distributions and fragmentation functions rise steeply towards small argument, this 
generally increases the relevance of the threshold regime and the soft-gluon effects are relevant 
even for situations where the the hadronic center-of-mass energy is much larger than the transverse 
momentum of the final state hadrons. This effect, valid in general in hadronic collisions, is even 
enhanced in single-inclusive hadron production since only a fraction z of the available energy is 
actually used to produce the final-state hadron. 



3 Resummed cross section 



We will now present the formulas for the threshold-resummed partonic cross sections. We will 
do this only for the case of the fully rapidity-integrated cross section, which turns out to signifi- 
cantly simplify the analysis. The resummation for the 77 dependence of the kinematically related 
prompt-photon cross section was performed in Ref. J7j, and we could in principle follow the 
techniques developed there to derive the 77 dependence of the resummed inclusive-hadron cross 
section. However, this process is of much greater complexity than prompt photons, as will become 
evident below, and it appears that a successful resummation at fixed rapidity will require further 
new techniques. We hope to address this in a future publication. As far as phenomenology is 
concerned, we will later on mimic the effects of the experimentally covered limited rapidity ranges 
by rescaling our resummed prediction by an appropriate ratio of NLO cross sections. Such an 
approximation was shown in J7] to work extremely well for the resummed prompt-photon cross 
section, where it was found that the shape of the cross section as a function of rapidity does not 
change much when going from NLO to the resummed result. This gives confidence that it may 
be applicable also in the case of inclusive-hadron production we are interested in. 



3.1 Mellin moments and threshold region 

The resummation of the soft-gluon contributions is carried out in Mellin- moment space, where 
the convolutions in Eq. (j2I) between parton distributions, fragmentation functions, and subprocess 
cross sections factorize into ordinary products. We take Mellin moments in the scaling variable 
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In A^-space Eq.© becomes 

a{N) = J2 fa/H, {N + 1, filj) h,H, {N + 1, ^il,) D,/e(2iV + 3, ^iIf) aab^,x{N) , (11) 



with the Melhn moments of the parton distribution functions and fragmentation functions, and 
where 
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Here, the threshold hmit w ^ 1 (or, for the rapidity-integrated cross section, x"^ —>■ 1) corresponds 
to — > oo, and the leading soft-gluon corrections arise as terms oc ag In^'^ A^. 

It is instructive to examine the interplay of rapidity integration and large- limit, for example 
in case of the NLO cross section in Eq. (jHJ. As we mentioned earlier, the soft-gluon terms are 
associated with the + distribution pieces in (jH}, which have coefficients that may be written as 
functions of v only. One has 



dv [4v{l - v)wf^^ f{v) 
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which implies that at large A^ the variable v is "squeezed" to v = 1/2, and hence, as follows 
from (j3]), the partonic rapidity is forced to f/ = 0. This means, for example, that near threshold 
it is justified to take the NLO partonic cross sections as proportional to the Born cross section. 
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with coefficients A',S',C" evaluated at v = 1/2. We will follow this reasoning also for the re- 
summed cross section to which we turn now. 



3.2 Resummation to NLL 

In Mellin-moment space, threshold resummation results in exponentiation of the soft-gluon cor- 
rections. Foremost, there are radiative factors for the initial and final partons, which contain 
the leading logarithms. At variance with the color-singlet cases of Drell-Yan and Higgs produc- 
tion 1^ 1^. and with prompt-photon production [T^ [T3] which has only one color structure 
at Born level, several color channels contribute to each of the 2 — > 2 QCD subprocesses relevant 
for inclusive-hadron production. As a result, there are color interferences and correlations in large- 
angle soft-gluon emission at NLL, and the resummed cross section for each subprocess becomes a 
sum of exponentials, rather than a single one. 
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In determining the resummed formula, we are in the fortunate situation that the effects of the 
color interferences for soft-gluon emission in the 2 — > 2 processes ab — *• cd have been worked out in 
detail in Refs. ^HlCSl for the case of jet production in hadronic collisions, which proceeds through 
the same Born processes. We take advantage of the formulas derived there. The results in [T8|[T9] 
have actually been given for arbitrary rapidity; for the case of the rapidity-integrated cross section 
we consider here it is sufficient to set r) = in the expressions of fHl ^| when diagonalizing 
(by changing the color basis) the soft anomalous dimension matrix computed there. A difference 
between inclusive hadrons and jets occurs regarding the treatment of the final-state parton c 
producing the jet or the hadron. In our case, this parton is "observed" , that is, we are considering 
a single-inclusive parton cross section. Such a cross section has final-state coUinear singularities 
which are factorized into the fragmentation functions. As far as resummation is concerned, the 
final-state observed parton therefore is similar to the initial-state partons and receives essentially 
the same radiative factor as the latter 



Combining these results of [TTl UHl HH 122] , we can cast the resummed partonic cross section 
for each subprocess into the rather simple form 



N '^N 
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^T-::liN) , (15) 



where the sum runs over all possible color configurations J, with G^^.^^^ representing a weight for 

each color configuration, such that X]/ ^ib-^cd — 1- ^'^Zcli-^) denotes the A^- moment expression 
for the Born cross section for the process, as defined in Eq. ()12j). We list the moment space 
expressions for all the Born cross sections in the Appendix. Each of the functions A^, J^, 
^{ir^)ab^cd -g^^ (ji5p is an exponential. represents the effects of soft-gluon radiation collinear 
to initial parton a and is given, in the MS scheme, by 

lnA%= / 4A>s(g^)), (16) 

and similarly for A^. Collinear soft-gluon radiation to parton c yields the same function, but 
with the initial-state factorization scale fipj replaced with the final-state one, Hff- The function 
Jff embodies collinear, soft or hard, emission by the non-observed parton d and reads: 
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/ / A^AMs{q')) + Mo^sil-z)Q') . (17) 

Jo 1-^ Ly(i-_^)2Q2 2 J 

Large-angle soft-gluon emission is accounted for by the factors ^(.'^)'^''^^'^^ which depend on the 
color configuration / of the participating partons. Each of the jg given as 

In Ai^?" = / -D,ab^cd{cis{{l - zfQ')) . (18) 

Jo i — 2 

Finally, the coefficient Cabled contains A^— independent hard contributions arising from one-loop 
virtual corrections. 

In the above formulas, Eqs. ()16 |) - (fTHj) . we have defined = 2p'^. Furthermore, each of the 
functions = Aa, B^, Dia.b->cd is a perturbative series in as, 

^(as) = ^-F«+f^)V(^) + ... , (19) 
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with [22]: 
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where Nf is the number of flavors, and 

Cg = CA = N, = ?> , C, = Cf = {Nl - l)/2iV, = 4/3 



7g = -3Cir/2 = -2 , 7g = -27r6o , &o 
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The expansion of the coefficients Cabled reads 

*cd 



C..^cA = 1 + ^C^Ld + 0{al] 



(21) 



(22) 



In the exponents, the large logarithms in now occur only as single logarithms, of the form 



Qfg In'^'*'^ ( A^) for the leading terms. Subleading terms are down by one or more powers of In(A^). 



Knowledge of the coefficients Aa 
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allows to resum the full towers of leading loga- 



rithms (LL) ag ln'^"'"^(A^), and NIL agln''(A^) in the exponent. Along with the coefficients C, 
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one then gains control of three towers of logarithms in the cross section, a^ln^^N) 



In 
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ab—>cd 

\N), 

(N), which is likely to lead to a much improved theoretical prediction. We also note that 
the factors A^y depend on the initial- or final-state factorization scales in such a way that they 
will compensate the scale dependence (evolution) of the parton distribution and fragmentation 
functions. One therefore expects a decrease in scale dependence, which indeed has been found in 
previous studies for other threshold-resummed cross sections. 

We finally examine the qualitative impact of the resummation. To this end, we note that, 
neglecting the running of the strong coupling, the LL terms in the exponents in Eqs. ()16|) and 
(fTTjl become 
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Therefore, for each partonic channel, the leading logarithms are 



^abScdW exp 
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The fact that this exponent is clearly positive for each of the partonic channels means that the soft- 
gluon effects will lead to an enhancement of the cross section. Particularly strong enhancements 
are to be expected for gluonic channels; for example, for the process gg — > gg one has Ca+Ch+Cc — 
Crf/2 = 15/2. The feature that partonic cross sections can give Sudakov enhancements is related 
to the fact that finite partonic cross sections are obtained after coUinear (mass) factorization, 
so that soft-gluon effects are partly already contained in the (MS-defined) parton distribution 
functions and, in our case, fragmentation functions. 



3.3 Exponents at NLL 



We now give explicit formulas for the expansions of the resummed exponents to NLL accuracy. 
Since the functions A^y and Jfj are "universal" in the sense that they depend only on the type of 
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the external parton, but not on the subprocess, their expansions are known and we recall them 
for the sake of completeness: 

lnA^(as(/^^),QV/^^;QV/^F) = ln7V/i«(A)+/ii^)(A,g7/x^;g7/.|,) + 0(as(«slniV)'=) , (25) 
lnJ^Mf^l),QVf^l) = ln7V/«(A) + /f (A,Q7//^) + 0(as(«sln7V)'=) , (26) 

where A = 6o'^s(A*/t) The functions h^^'"^^ and f^^'"^^ are given by 
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Here, as before bo = (HCa - 2Nf) /Un, and 
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In ^ . 
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We remind the reader that the scale fip represents the initial-state (final-state) factorization scale 
IJ^Fi {i^ff) for the radiative factors for the initial (final) state. The functions h^^^ and f^^^ above 
contain all LL terms in the perturbative series, while /i^^) ^nd /*^2) qj-q Qf nll accuracy only. For 
a complete NLL resummation one also needs the coefficients In /^(.^^)"''~*"^'^ whose NLL expansion 
reads: 



In A 
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VAT 

As we mentioned earlier, the D^'lh^^.^i the corresponding "color weights" Gjab~*cd, ^.tc both 



process and "color configuration" dependent. All the coefficients -Oj^j^^^ and Giab^cd that we 
need to NLL are listed in the Appendix. 
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3.4 Coefficients C);^Ucd 

We have verified for eacli subprocess tliat expansion of tlie resummed formulas above to 0{q^) 
correctly reproduces the logarithmic terms oc ag In^(A^), ag In(iV) known from the full fixed-order 
calculations |21 E] • Comparison to those calculations also allows to extract the first-order coeffi- 
cients C^ah^cd- Numerical results for the coefficients are presented in the Appendix. 



3.5 Matching to the NLO cross section, and inverse MelUn transform 

As we have discussed above, the resummation is achieved in Mellin moment space. In order to 
obtain a resummed cross section in x\ space, one needs an inverse Mellin transform. This requires 
a prescription for dealing with the singularity in the perturbative strong coupling constant in 
Eqs. fll6|) - (fTSj) or in the NLL expansions, Eqs. (j7fj) - (jH^ . We will use the Minimal Prescription 
developed in Ref. [21], which relies on use of the NLL expanded forms Eqs. (jTTjl - ljH^ . and on 
choosing a Mellin contour in complex- space that lies to the left of the poles at A = 1/2 and 
A = 1 in the Mellin integrand: 

= f""- ^JL (4)-".(».(iV) , (33) 

dpT JcMP-ioo 27r2 

where 6otts(/^|;) IhCa/p < 1/2, but all other poles in the integrand are as usual to the left of the 
contour. The result defined by the minimal prescription has the property that its perturbative 
expansion is an asymptotic series that has no factorial divergence and therefore no "built-in" 
power-like ambiguities. Power corrections may then be added, as phenomenologically required. 

When performing the resummation, one of course wants to make full use of the available 
fixed-order cross section, which in our case is NLO {0{a\)). Therefore, a matching to this cross 
section is appropriate, which may be achieved by expanding the resummed cross section to 0(Q;g), 
subtracting the expanded result from the resummed one, and adding the full NLO cross section: 



^3 ^^(match)^^^) rCMP+ioo ^.-N+l 



dpi 



OK) 



dpT 



where cr'^J^b^cdi-^) resummed cross section for the partonic channel ab — > cd as given in 

Eq. p5|) . In this way, NLO is taken into account in full, and the soft-gluon contributions beyond 
NLO are resummed to NLL. Any double-counting of perturbative orders is avoided. 



4 Phenomenological Results 

Starting from Eq. (jS^), we are now ready to present some first resummed results at the hadronic 
level. This is not meant to be an exhaustive study of the available data for inclusive-hadron 
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production; rather we should hke to investigate the overall size and relevance of the resummation 
effects. We will only consider 7r° production in pp collisions and compare to a few selected sets of 
data. 

Let us begin by specifying some "default" choices for the distribution functions that we will 
use in our studies. We will use the MRST2002 set of parton densities and the pion frag- 
mentation functions of j2E] (referred to as "KKP"). For comparison, we will also present some 
results for Kretzer's set j^Z] of fragmentation functions. Note that, according to Eq. ()34|) . it is a 
great advantage to have parton densities and fragmentation functions available in moment space. 
Technically, since the MRST distributions are not available in moment space, we first performed 
a fit of a simple functional form to the MRST distributions, of which we were then able to take 
moments. This had to be done separately for each parton type and for each scale. Concerning the 
fragmentation functions, Kretzer's set is anyway set up in moment space, and we found it possible 
to analytically take moments of the KKP parameterization. 

As we discussed at the end of subsection 13.21 we generally expect fairly large effects from 
soft-gluon resummation for inclusive-hadron production. This makes it rather important to be 
sure that the resummed soft-gluon terms indeed constitute the dominant part of the cross section 
and do not, for example, lead to an overestimate of the higher orders. We therefore start by 
identifying the kinematic regions where soft-gluon contributions are likely to dominate the cross 
section. A gauge for this is obtained by comparing the resummed formula expanded to NLO to the 
full fixed-order (NLO) perturbative result, that is, by comparing the last two terms in Eq. ()34|) . 
Figure HI shows this comparison for a typical fixed-target energy = 31.5 GeV, and for RHIC's 
^/S = 200 GeV. As can be observed, the expansion faithfully reproduces the NLO result. In 
the fixed-target regime the agreement is excellent, except perhaps at the lowest pion transverse 
momenta, ~ 3 GeV, where the soft approximation tends to yield a slight overestimate. This 
is obviously related to the fact that the smaller px (at fixed energy), the further one is away 
from threshold, so that the soft-gluon approximations become less reliable. The same is expected 
to happen if the energy is increased at fixed transverse momentum. Indeed, as the curves for 
■\/s = 200 GeV in Fig. Q show, the NLO-expanded resummed result, while still remarkably close 
to the full NLO prediction, gives a less accurate picture of the latter than at fixed-target energies. 
Our conclusion from Fig. ^ is therefore that the contributions associated with the near-threshold 
region are dominant in the fixed-target regime, implying that resummation will be relevant and 
accurate here, even at relatively small transverse momenta. At colliders, our resummed cross 
section will likely be too large, and further improvements in the theoretical framework may be 
needed^. We will briefly return to this point at the end of this paper. In the following we will 
primarily focus on the fixed-target regime. 

We next investigate how large the higher-order contributions provided by NLL resummation 
are. To this end, we go back to Eq. (jH^ and take the full resummed result, defined in the minimal 
prescription and matched to NLO. As before, the cross sections are integrated over all rapidities. 
We define a resummed "i^-factor" as the ratio of the resummed cross section to the NLO cross 
section, 

~ da^^^^)/dpT ' ^ ^ 

tit is worth recalling that our matching procedure given by Eq. (|34|) ensures that the NLO cross section is 
always fully and exactly taken into account in our final "matched" cross section, so that any overestimate would 
only occur at NNLO and beyond. 
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Figure 1: Comparison of full NLO cross sections (solid lines) for pp vr^X with the NLO {0{a^)) 
expansion of the resummed cross section (dashed lines), for two different energies. We have chosen 
the factorization and renormalization scales as pt- 



which is shown for the fixed-target regime, and for scales hr = hfi = IJ-ff = Pt, by the solid line 
in Fig. 1^ As can be seen, 

^(res) 

is very large, meaning that resummation results in a dramatic 
enhancement over NLO. It is then interesting to see how this enhancement builds up order by 
order in the resummed cross section. We therefore expand the matched resummed formula beyond 
NLO and define the "soft-gluon i^-f actors" 

= dai^^o)/dpT ' ' ^^^^ 

which for n = 2, 3, . . . give the additional enhancement over full NLO due to the 0{a^^^) terms in 
the resummed formula. Formally, = 1 and K°° = K^^^^^ of Eq. (jHSj) . The results for 7^2,3,4,5,6 
are also shown in Fig. |21 One can see that there are very large large contributions even beyond 
NNLO, in particular at the higher px- Clearly, the full resummation given by the solid line is 
required here. 

As we have mentioned earlier, we have determined the resummed formulas for the fully rapidity- 
integrated cross section, whereas in experiment always only a certain limited range of rapidity is 
covered. In order to be able to compare to data, we therefore approximate the cross section in 
the experimentally accessible rapidity region by 

p3^^{match) _ 3^^^(NLO) _ 

; (77 mexp. range) = il ' ; (r? mexp. range) , (37) 

dpT dpT 

where K^"^*^^^ is as defined in Eq. (I35|) in terms of cross sections integrated over the full region 
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Figure 2: "A'-factors" relative to NLO as defined in Eqs. ()35|) and ()36p for pp n^X in tfie 
fixed-target regime. 



of rapidity. In otlier words, we rescale tfie matcfied resummed result by the ratio of NLO cross 
sections integrated over the experimentally relevant rapidity region or over all rj, respectively. 

In Fig. El we compare the NLL resummed and NLO predictions to the available data from 
E706 I2H] for pp -> TT^X at VS = 31.5 GeV. The data cover |r/| < 0.75. We use the KKP 
fragmentation functions j2ni and give results for three different choices of scales, fiR = Hfi = 
l^FF = CPt, where C = 1/2, 1, 2. It is evident that the NLO result falls far short of the data, which 
is an observation that has been made before El EI- Furthermore, there is a very large scale 
dependence at NLO. The situation is significantly improved when the NLL resummation is taken 
into account. As we already saw in Fig. the NLL matched cross section is considerably higher 
than the NLO one, and it shows a markedly improved comparison to the E706 data, probably 
satisfactory in view of the overall uncertainties. Furthermore, the scale dependence is considerably 
reduced compared to the NLO calculation, and hence the accuracy of the prediction is improved. 

FigureEJshows the same result, but now for the Kretzer set |2Z] of pion fragmentation functions. 
These functions are known to be overall significantly smaller than the ones of KKP, in particular 
for the gluon fragmentation function which is not well determined from the e~^e~ hX data^. 
One therefore finds that all theory curves are shifted downward with respect to the results shown 

■l-Note that analyses of hadron production in the additional jet in e+e^ — ^ 66jet events (2S! do constrain Dg 
significantly. The Dg 's in the sets of and |22j are in reasonable agreement with these data, with the one of P?) 
arguably setting a lower bound on Dg. 
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Figure 3: NLO and NLL resummed results for the cross section for pp —>■ n^X for E706 kinematics. 
We have used the KKP fragmentation functions I^Hl- Results are given for three different choices 
of scales, fiR = fipi = fJ'FF = CPt, where ( = 1/2, 1, 2. Data are from ^Hj . 



in the previous figure. Nevertheless, the effects due to threshold resummation remain large. 

To give another example. Fig. El presents a comparison to data from the WA70 experiment, 
corresponding to = 31.5 GeV and \xp\ < 0.45. Again a significant enhancement due to NLL 
resummation is found, resulting in a much improved agreement between theory and data. 

Finally, in Fig. IHl we repeat the calculations for the case of proton-proton collisions at RHIC 
with = 200 GeV and \ri\ < 0.35. The data are from the measurement performed by the 
PHENIX Collaboration [Hj. Again, an enhancement from resummation is found which however is 
smaller than in the previous figures. This is expected since we are further away from threshold 
here, due to the much higher energy. Nevertheless, the enhancement is quite significant at the 
larger pt, where in fact the resummed result appears to lie too high. We emphasize, however, that 
according to our results shown in Fig.^it is likely that the NLL resummation gives an overestimate 
of the higher-order corrections in this case. We therefore do not take the enhancement too literally 
and reserve its closer investigation to a future study. We note that also in this case there is a 
considerable reduction in the scale dependence, and that again the fragmentation functions of [2Z] 
lead to a smaller cross section. 
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Figure 4: Same as Fig. |21 but for the set of fragmentation functions of j2Z] ■ 



5 Conclusions and outlook 

We have studied the NLL all-order resummation of threshold logarithms in the partonic cross 
sections relevant for the process pp —>■ hX at high transverse momentum of the hadron h. This 
study has in part been motivated by the observed shortfall of fixed-order (NLO) cross sections when 
compared to experimental data in the fixed-target regime, in contrast with the excellent agreement 
of data and theory at colliders. Our numerical results indeed show a strong enhancement of the 
cross section over the next-to-leading order one for typical fixed-target kinematics, significantly 
improving the agreement between data and theoretical predictions. At higher energies, such as at 
RHIC, the resummation effects are less important, but more theoretical analysis is needed here 
due to the likely relevance of subleading terms. 

We emphasize that the contributions generated by resummation are a well-defined class of 
higher-order corrections to the leading-power partonic cross sections that will be present in the 
full perturbative series order by order and actually dominate it. Our results are then also to be seen 
in the context of the size of possible non-perturbative power-suppressed corrections to the cross 
section. Any residual shortfall of the resummed theoretical prediction would need to be attributed 
to such contributions. In previous studies, "intrinsic" transverse momenta of partons have often 
been taken into account in (LO or NLO) calculations of inclusive-hadron cross sections jHlE]; in 
order to bridge the large gaps between data and NLO theory in the fixed-target regime. These 
can perhaps best be viewed as models of the power-suppressed contributions. In the light of our 
results, however, much of the enhancements needed for a satisfactory description of the fixed-target 
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Figure 5: Same as Fig. El but comparing to the pp n^X data of WA70 [^0] at ^/S = 22.9 GeV. 



data appears to come from perturbative contributions, so that power-suppressed contributions are 
probably of rather moderate size. It is interesting to note that resummed perturbation theory 
itself may provide information on the structure of power corrections, through ambiguities in the 
perturbative series arising from the pole in the perturbative running coupling in the expressions 
Eqs. flT^ - (fTHj) . A recent study addressed this issue in the case of single-inclusive cross sections 
at large xt and indeed estimated power corrections to be not very sizable. On the other hand, it 
is known IT?)1 ITT] that threshold resummation effects are not very large for the prompt photon 
cross section in the fixed-target regime, where discrepancies between data and NLO theory of 
similar magnitude as for pion production have been observed in some cases. This issue clearly 
needs further study. 

We finally emphasize that we regard this study only as the beginning of a more detailed anal- 
ysis of threshold resummation for inclusive-hadron cross sections. There are several points in 
which further developments are desirable. First of all, as we noted earlier, it would be possible 
in principle (albeit challenging technically) to perform the resummation correctly for the fully 
rapidity-dependent cross section. To do this appears all the more interesting since it was ob- 
served that the discrepancies between NLO and fixed-target data actually increase (at fixed 
Pt) toward larger rapidities. We expect resummation effects to become even more important as 
well at large t], simply because one is approaching threshold more closely. It is also possible to 
improve the resummation by resumming also terms of the form a^\n'^^~^{N)/N in the partonic 
cross sections. Such terms arise from coUinear emissions [201 1^ 121] ! they are suppressed with 
respect to the LL and NLL terms but may nonetheless be of relevance if one is further away 
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Figure 6: Same as Fig. El but comparing to the RHIC PHENIX data of p. 



from threshold as for example in the fixed-target regime at lower p^, or a collider energies. We 
expect that such additional contributions would also further decrease the scale dependence. With 
these improvements in place, detailed phenomenological studies of power corrections might be- 
come possible. We finally also note that another important field for further study would be the 
effects of threshold resummation on spin asymmetries, in particular on the double-longitudinal 
spin asymmetries An measured by E704 [33] and at RHIC |l9|. We beheve that the significant 
enhancements due to resummation that we found in this work strongly motivate all these further 
studies. 
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Appendix: Results for the various subprocesses 



In this appendix we compile the moment-space expressions for the Born cross sections for the 
various partonic subprocesses, and the process dependent coefficients cf),!^^^, Diab^cd, and G^^^^^. 
Since the C^^^^d have rather lengthy expressions, we only give their numerical values for Nf = 5 
and the factorization and renormalization scales set to fj^pj = fj^pp = fJ^R = Q- 



qq' qq': 



.(Born) 
qq'^qq 



.(N) ^^{5N^ + 15N + 12) b(^N,1^ ' 



qq' qq' 



Glqq'^qq' - 1/3 , G2qq'^qq' - 2/3 , D^lgq'-^gq' - "4 lu 2 , D^2qq'-^qq' " ^ ' 

^il'-^,,' = 20.2389 {Nf = 5) . (38) 



-^q'i^) - 1^ (5iV^ + 15iV + 12) B (iV, , 

G,qq^-.qq>-l/9, = 8/9 , i^^ =-10/3 1n2 D^;i-, -,^8/3 \n2, 



lqq—*qq ' 2qq'—*qq' 

C;;)^-,^^^-, = 22.4483 {Nf = 5) . (39) 



qq q'q': 



Glqq^q'q' = 1 , /^j^J-^^,-, = -10/3 lu 2 , , 

C'f.U'.' = ^-91881 {Nf = 5) . (40) 



qq qq: 



qq—^qq \ J 



qq^qqy I 'iC\ 



(Ca(5A^2 + 15A^ + 12) - 2A^(3 + 2N)) B [n , , 



Glqq^qq — 9/11 , G2qq-^qq — 2/11 , D^\lq-^qq — "4 lu 2 , D^lq^qq — , 

Cf^U. = 19-535 {Nf = 5) . (41) 



qq — > qq: 



a^SZi{N) = ^ {Ca{IIN' + 59iV2 + 102iV + 60) + N{N + 3)(5 + 2N)) B (^iV, , 

Gi,,-_,,- = 5/21, G2M = 16/21, i5S,-_,,-=-10/3 1n2, L>g,-^^^- = 8/3 In 2 , 
C?qq-.qq = 19.9643 (7V^ = 5) . (42) 
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m 99- 



Q9 Q9- 



qg gq: 



99 99' 



99 QQ- 



a'^Z7o{N) = ^ {2Cf{N + 2)(5 + 2N) - Ca{N +1){N + 3)) Bi^N + 1,- 
(71,,^,, = 5/7, G'2,,_,, = 2/7, <),^,, = -10/3 1n2, L>«-_^^ = 8/3 In 2 , 



M = 12.4329 {Nf = 5) . 



^^J^Z^giN) = ^ {CfN{7 + 57V) + 2C^(57V2 + 157V + 12)) B (^N, ^) 

Glqg^qg — 45/88 , G2qg^qg — 25/88 , G^qg-^qg — 18/88 , 

D^'lg^^g = -14/3 In 2 , D'il^^^ = 10/3 In 2 , d'H^^^ = -2/3 In 2 , 
Cf^,^,, = 15.4167 {Nj = 5) . 



aJJ^^°J(A^) = ^ (Cf 7V(7 + 57V) + 2Ca{5N^ + 157V + 12)) B (^N, 

Glqg^gq = 45/88 , G2qg^gq = 25/88 , Gsgg^gq = 18/88 , 

<U. = -8 1^2, <U, = 0, <,^,, = -4 1n2, 
Cil^gg = 22.4474 (7V^ = 5) . 



^SZiiN) = ^ (97V' + 457V2 + 727V + 40) B (n, l) 

G^gg^gg = V3 , G2gg^gg = 1/2 , G^gg^gg = 1/6 , 

ci;U. - 21.1977 (7V, = 5) . 



= (2C^(A^ + 2)(5 + 27V) - C47V + 1)(7V + 3)) ^ (^7V + 1, ^ 

G^l gg~*qq = 5/7 , G2gg^qq = 2/7 , D^lgg-,qq = , D^^^lg^qq = 6 In 2 , 

M = 16.7962 (7V^ = 5) . 
In the above expressions, B{a, h) is the Beta-function. 
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